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INTRODUCTION 


A  hollow  cylinder  Is  a  common  configuration  used  in  pressure  vessels. 
According  to  the  classical  Lame"  solution  for  an  elastic  cylinder,  the  largest 
tangential  stress  occurs  at  the  bore.  In  order  to  increase  the  maximum  pres¬ 
sure  a  cylinder  can  contain,  to  reduce  the  probability  rate  of  crack  initia¬ 
tion  at  the  bore  and  to  slow  down  the  growth  rate  of  a  fatigue  crack,  a  resid¬ 
ual  compressive  stress  can  be  introduced  near  the  internal  boundary.  ”Auto- 
frettage”  is  such  a  process  used  to  produce  the  favorable  residual  stresses  in 
a  gun  barrel.  By  swaging  an  oversized  mandrel  through  the  bore  or  by  apply¬ 
ing  a  sufficiently  high  hydrostatic  pressure  at  the  bore,  the  cylinder  is  par¬ 
tially  yielded.  Upon  removal  of  the  inner  pressure,  a  residual  stress  distrl^ 
bution  will  result  in  the  cylinder.  The  tangential  stress  is  highly  compres¬ 
sive  at  the  bore  and  varies  logarithmically  to  tension  over  the  region  which 
was  plastically  deformed  during  pressurization.  The  plane-strain  residual 
stresses  due  to  autofrettage  are  well-known.  The  residual  stress  pattern  will 
change  due  to  changes  of  geometrical  configurations  such  as  the  presence  of 
keyways,  notches,  cracks,  etc.  The  redistribution  of  residual  stresses  in 
these  Instances  is  not  generally  straight  forward. 

In  a  previous  article^  we  showed  that  an  active  thermal  load  can  be  used 
to  produce  thermal  stresses  equivalent  to  autofrettage  residual  stresses. 

This  approach  provides  a  method  to  find  the  redistribution  of  residual 


^Hussain,  M.  A.,  Pu,  S.  L. ,  Vasllakis,  J.  D. ,  and  O'Hara,  P. ,  ’’Simulation  of 
Partial  Autofrettage  by  Thermal  Loads,”  Journal  of  Pressure  Vessel 
Technology,  Vol.  102,  No.  3,  1980,  pp.  314-318. 


stresses  by  computing  the  thermal  stresses  In  the  cylinder  with  notches  and 
cracks  under  the  thermal  load.  In  this  report  we  supply  numerical  verifica¬ 
tion  to  show  the  equivalence  of  thermal  and  residual  stress  redistribution 
caused  by  notches  and  the  equivalence  in  stress  intensity  factors  for  cracks 
due  to  residual  stresses  and  simulated  thermal  loads.  The  method  of  super¬ 
position  requires  distributed  loads  on  the  crack  face.  It  results  in  singular 
stresses  near  the  crack  tip  from  the  finite  element  computations  using  APES. 

On  the  other  hand,  the  method  of  thermal  simulation  is  simple  in  application 
and  yields  accurate  results  including  stresses  near  the  crack  tip. 


THERMAL  STRESS  SIMULATION 

The  plane-strain  stress  distribution  of  a  partially  autofrettaged  tube 
using  the  von  Mises’  yield  condition  and  the  incompressibility  condition  is 


given  by2 


°o  r  Pz  l  l  . 

—  {(2  log  -  -  1  +  -r)  -  Pi (-j  -  -5)}  a  <  r  <  p  (1) 
/3  P  bz  bz  rz 


~  <P2-Pl)(-2  -  -5> 


p  <  r  <  b  (2) 


°o  .  r  p2  11, 

—  {2  log  -  +  1  +  -=■  -  Pi (-5-  +  -5)}  a  <  r  <  p  (3) 

/3  p  bz  bz  rz 


°o  ,  11 

n  ““  Pl)(P  +  ?’ 


P  <  r  <  b  (4) 


2Hill,  R..  The  Mathematical  Theory  of  Plasticity.  Oxford  at  the  Clarendon 
Press,  1950. 
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Where  P},  a  constant  depending  on  geometrical  dimensions  only 


P1 


a2b2 

b2-a2 


Cl 


p‘  p 

?  +  2  108  ;> 


(5) 


Upon  the  application  of  the  following  thermal  load 

(Ta-Tp) 
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the  thermal  stresses  In  the  hollow  cylinder  are* 
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^Hussain,  M.  A.,  Pu.  S.  L. ,  Vasllakls,  J.  D. ,  and  O'Hara,  P.,  "Simulation  of 
Partial  Autofrettage  by  Thermal  Loads,”  Journal  of  Pressure  Vessel 
Technology,  Vol,  102,  No.  3,  1980,  pp.  314-318. 
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If  Che  temperature  gradient  of  the  thermal  load  and  the  yield  stress  of  the 
cylinder  material  satisfy  the  relation 


Ea(Ta-Tp)  2 0o 

2(l-v)log(p/a)  /3 


(8) 


the  thermal  stresses  and  the  autofrettage  residual  stresses  become 
equivalent . 

Using  the  temperature  distributions  as  temperature  Input  In  a  finite-dif¬ 
ference  computer  program  based  on  the  theory  of  thermal  stress  in  Sections  9- 
10  of  reference  3,  numerical  results  were  obtained  for  the  thermal  stresses  in 
a  cylinder.  These  values  were  shown  in  reference  1  to  be  in  good  agreement 
with  residual  stresses  computed  from  equations  (1)  to  (4).  Thermal  stresses 
were  also  computed  using  the  finite  element  NASTRAN  program.  The  finite  ele¬ 
ment  results  for  the  thermal  stresses  were  also  In  excellent  agreement  with 
the  autofrettage  residual  stresses. 1  It  remains  to  be  shown  that  the  redis¬ 
tribution  of  thermal  stresses  due  to  the  presence  of  geometrical  discontinu¬ 
ities  simulates  the  redistribution  of  residual  stresses  caused  by  the  same 
geometrical  changes.  This  question  of  simulation  was  specifically  raised  by  a 
review  of  the  work  of  Kapp,^  who  applied  the  idea  of  thermal  load  simulation 
to  stress  analysis  of  autofrettaged  cylinders  with  an  0D  notch. 


^Hussain,  M.  A.,  Pu,  S.  L. ,  Vasilakis,  J.  D. ,  and  O'Hara,  P.,  "Simulation  of 
Partial  Autofrettage  by  Thermal  Loads,"  Journal  of  Pressure  Vessel 
Technology,  Vol.  102,  No.  3,  1980,  pp.  314-318. 

^Boley,  B,  A.  and  Weiner,  J.  H. ,  Theory  of  Thermal  Stresses,  John  Wiley  and 
Sons,  1960. 

^Kapp,  J.  A.  and  Pflegl,  G.  A.,  "Stress  Analysis  of  0D  Notched  Thick-Walled 
Cylinders  Subjected  to  Internal  Pressure  or  Thermal  Loads,”  Technical  Report 
ARLCB-TR-80005,  Benet  Weapons  Laboratory,  LCWSL,  ARRADCOM,  US  Army,  February 
1980. 
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STRESS  INTENSITY  FACTOR  DUE  TO  THE  RESIDUAL  STRESS 

Consider  a  crack  in  Che  radial  direction  emanating  from  the  outer  bound¬ 
ary  of  an  autofrettaged  tube.  The  crack  will  be  opened  up  due  to  the  tensile 
residual  stress  in  the  circumferential  direction.  This  will  change  the  resid¬ 
ual  stress  distribution  in  the  tube.  A  method  is  required  to  compute  the 
stress  intensity  factor  due  to  the  residual  stress.  The  method  of  thermal 
simulation  was  proposed  and  used  in  reference  1.  The  method  of  superposition, 
which  will  be  discussed  below,  is  another.  The  result  from  the  method  of 
superposition  is  within  one  percent  of  that  reported  in  reference  l  by  the 
method  of  thermal  simulation. 

For  an  ID  (inner  diameter)  radial  crack,  the  crack  remains  closed  in  the 
compressive  residual  stress  region.  The  stress  intensity  factor  is  zero  and 
the  crack  does  not  affect  the  residual  stress  field.  The  crack  will  open  when 
a  sufficiently  large  internal  pressure  is  applied  on  the  inner  boundary  of  the 
cylinder.  The  stress  intensity  factor  is  the  net  sum  of  two  values,  one  is 
positive  corresponding  to  the  tensile  stress  due  to  the  applied  internal 
pressure  and  the  other  is  negative  due  to  the  compressive  residual  stress.  If 
the  absolute  value  of  the  latter  is  larger  than  the  former,  the  net  result  of 
the  stress  intensity  factor  is  zero  since  negative  stress  intensity  factor  has 
no  physical  meaning.  However  for  computational  purposes,  a  negative  stress 
intensity  factor  due  to  the  compressive  residual  stress  is  convenient  and 
should  be  understood  as  such.  It  is  necessary  to  have  an  effective  method  for 

^Hussain,  M.  A.,  Pu,  S.  L. ,  Vasilakls,  J.  D. ,  and  O'Hara,  P. ,  "Simulation  of 
Partial  Autofrettage  by  Thermal  Loads,"  Journal  of  Pressure  Vessel 
Technology,  Vol.  102,  No.  3,  1980,  pp.  314-318. 
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the  computation  of  the  stress  intensity  factor  due  to  the  compressive  residual 
stress. 

The  thermal  stress  simulation  mentioned  in  the  previous  section  provides 
a  method  of  computation  of  stress  intensity  factors  due  to  residual  stresses. 
The  residual  stress  is  replaced  by  an  active  thermal  load  which  produces  a 
thermal  stress  distribution  identical  to  the  residual  stress.  The  stress 
intensity  factor  sought  is  obtained  from  that  due  to  the  thermal  load. 

The  finite  element  computer  program  APES,  using  12-node  isoparametric 
elements,  is  used  in  all  the  numerical  computations  reported  here.  This 
powerful  program  has  been  continually  expanded  to  add  new  features  such  as 
graphical  output,  thermal  loading,  crack  face  loading  for  enriched  elements, 
etc.,  and  to  increase  the  capacity  to  handle  a  maximum  of  200  elements. ^ 

The  new  thermal  input  feature  of  APES  came  just  in  time  for  the  fracture  anal¬ 
ysis  due  to  thermal  loading. 

As  a  first  example,  consider  a  cylinder  with  two  ID  cracks,  Figure  1(a). 
From  the  symmetry  of  the  problem,  only  the  first  quadrant  of  the  cylinder  is 
needed  to  be  considered.  The  finite  element  idealization  using  two  enriched 
elements  at  the  crack  tip  is  shown  in  Figure  1(b).  The  temperature  at  the 
bore  is  arbitrarily  chosen  to  be  Ta  •  0.  The  temperature  at  the  outer  bound¬ 
ary  can  be  computed  from  equation  (8).  It  depends  on  the  material  constants 

^Gifford,  L.  Nash,  "APES-Finite  Element  Fracture  Mechanics  Analysis  Revised 
Documentation,”  Naval  Ship  Research  and  Development  Center  Technical  Report 
DTNSRDC-79/023,  March  1979. 
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of  the  cylinder  as  well  as  the  degree  of  overstrain.  For  a  cylinder  with  E  ■ 
30  x  106  psl  (207  GPa) ,  v  «  0.3,  a  »  6.8  x  10”6  in/in/®F  (12.24  x  10”6  cm/cm/ 
®C),  0o  ■  170  Ksi  (1173  MPa),  the  temperature  at  the  outside  boundary  Is 
Tp  *  -633.17°F  (-333.98°C)  for  60%  autof rettage.  The  same  quantity  for  100% 
autofrettage  is  Tp  =  -933.77°F  (-500.98°C) .  For  a  crack  depth  d  *  0.2  in.  In 
a  tube  of  b  *  2  in.,  a  =  1  in.  the  numerical  results  for  the  stress  intensity 
factor  due  to  the  thermal  load  are  Kj  =  -89.403  Ksi/in  (1  Ksi  /Tn  =1.1  MPa^m) 
and  Ki  *  -109.63  Ksi/Tn  for  60%  and  100%  overstrained  tubes,  respectively. 
Using  the  well-known  form  of  Kj  *  l.  12(ird)  ^^0o(r=a)  for  shallow  edge  cracks 
as  a  normalizing  factor,  the  corresponding  dimensionless  stress  intensity 
factors  are  Kj/[ 1. 1 2 ( wd) */2(a0)a]  =  0.70  and  0.74,  respectively. 

Another  method  based  on  the  principle  of  superposition  can  be  used  to 
compute  the  stress  intensity  factor  due  to  the  autofrettage  residual  stress 
denoted  by  SR.  The  tangential  residual  stress  of  SR  on  the  cross  section  c-c. 
Figure  1(a),  for  a  100%  autofrettaged  cylinder  is  shown  in  Figure  2(a).  The 
stresses  at  r  ■  a,  a+d,  and  b  are  denoted  by  (oe)a,  (oe)d  and  (°0)b»  respec¬ 
tively.  Let  Sc  be  the  normal  traction  distributed  on  the  crack  surface  which 
tends  to  close  the  crack.  Its  magnitude  is  identical  to  the  corresponding 
tangential  residual  stress  of  SR  as  shown  in  Figure  2(b).  The  normal 
traction  which  is  equal  and  opposite  to  Sc  is  designated  by  S-p,  Figure  2(c). 

a 

If  both  Sc  and  Sx  are  added  to  SR,  they  produce  no  net  effect,  Figure  3(a). 

Now  if  Sx  and  SR  are  combined,  the  crack  surface  is  free  of  stress  and  remains 
closed.  The  stress  intensity  factor  is  of  course  zero,  Figure  3(b).  The  neg¬ 
ative  stress  intensity  factor  due  to  the  normal  traction  Sc ,  Figure  3(c),  must 


equal  the  stress  Intensity  factor  due  to  Sr.  The  same  finite  element  ideal¬ 
ization  of  Figure  1(b)  is  used  in  the  computer  program  to  obtain  the  stress 
intensity  factor  due  to  Sc,  which  replaces  the  thermal  load  in  the  previous 
computation.  The  new  feature  of  the  APES  program  of  crack  face  loading 
for  enriched  elements  is  used.  The  numerical  result  for  a  distributed  load 
Sc  corresponding  to  the  100%  overstrain  is  Kj  -  -109.44  Ksi/in.  The  percent¬ 
age  of  error  for  the  previous  result  due  to  thermal  loading  is  0.17%.  This 
confirms  that  the  thermal  simulation  gives  a  method  for  the  computation  of  the 
stress  Intensity  factor  due  to  residual  stress. 

The  APES  program  with  the  distributed  crack  face  loading  gives  the  cor¬ 
rect  stress  Intensity  factor  and  singular  stresses  at  the  crack  tip  and  its 
adjacent  nodes.  For  instance,  the  stress  og  at  the  node  nearest  to  and  in 
front  of  the  crack  tip  along  the  x-axis  is  -0.786  x  104^  psi.  The  thermal 
stress  at  the  same  node  due  to  the  simulated  thermal  load  is  og  -  -0.564  x  106 
psi.  From  the  stress  Intensity  factor  obtained,  the  normal  stress  can  be 
estimated  by  the  near  field  solution  due  to  Westergaard^ 

*1  0  0  36 

o,  ”  -  cos  -  [1  -  sin  -  sin  — ]  (9) 

✓7F  2  22 

K1  0  6  30 

ov  ■  -  cos  -  [1  +  sin  -  sin  — ] 

y  2  2  2  (10) 

°Westergaard ,  H.  M. ,  "Bearing  Pressures  and  Cracks,”  Transactions  of  the  ASHE, 
Journal  of  Applied  Mechanics,  1939. 
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For  9  *  0,  r  *  0.01S  and  ■  -109.44  Ksi/in,  we  have  09  ■  Oy  *  -0.504  x  106 

pal.  The  method  of  thermal  simulation  gives  not  only  correct  stress  Intensity 
factors  for  cracks,  but  also  correct  stresses  even  close  to  the  crack  tip. 

STRESS  REDISTRIBUTION 

Due  to  the  presence  of  notches,  such  as  the  two  notches  shown  In  Figure 
4(a),  redistribution  of  residual  stresses  occurs  in  the  cylinder.  To  find  the 
resultant  residual  stresses,  thermal  stress  simulation  may  be  used.  The  redis¬ 
tribution  of  thermal  stresses  due  to  the  presence- of  notches  under  the  simu¬ 
lated  thermal  load  should  be  the  same  as  residual  stress  redistribution.  This 
will  be  verified  numerically  in  what  follows. 

The  redistribution  of  thermal  stresses  can  be  obtained  numerically  by  the 
computer  program  APES.  The  finite  element  mesh  of  a  quarter  of  the  ring  with 
one  half  of  the  notch  is  shown  in  Figure  4(b).  The  tangential  thermal  stresses 
at  nodal  points  1-10  along  the  x-axis  and  at  nodal  points  88-100  along  the  y- 
axis  are  tabulated  in  column  (1)  of  Table  I  for  a  temperature  gradient  corre¬ 
sponding  to  60%  overstrain. 

The  redistribution  of  residual  stresses  can  be  computed  based  on  the  prin¬ 
ciple  of  superposition.  Let  Sg  and  Sg'  be  the  residual  stresses  in  the  ring 
before  and  after  the  introduction  of  the  notches.  The  compressive  tangential 
stress  on  mn  and  the  compressive  radial  stress  on  nq  can  be  computed  from  equa¬ 
tions  (3)  and  (1)  for  the  ring  without  the  notch.  The  state  of  stress  produced 
in  the  ring  due  to  the  distributed  compressive  stresses  just  obtained  on  mn  and 
nq  is  denoted  by  Sc.  Let  be  equal  and  opposite  to  Sc.  The  redistribution 
of  residual  stresses,  Sg' ,  is  not  affected  by  adding  Sc  and  S?  simultaneously. 
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The  sum  of  Sr'  and  Sc  Is  nothing  but  Sr.  Hence  Sr*  is  the  sum  of  Sr  and  Sf. 
Exact  values  of  tangential  residual  stress  for  Sr  at  the  selected  nodal  points 
are  tabulated  in  column  (2)  of  Table  I.  The  state  of  stress  S^  can  be  com¬ 
puted  by  APES  using  the  same  finite  element  mesh,  Figure  4(b),  with  the 
applied  load  distributed  on  mn  and  nq.  The  numerical  results  of  a q  for  Sf  are 
given  in  column  (3)  of  Table  I.  The  sum  of  two  values  in  columns  (2)  and  (3) 
is  eg  for  Sr'  given  in  column  (4).  The  agreement  between  09  in  column  (4)  for 
Sr'  and  that  in  column  (1)  for  thermal  load  is  excellent. 

Table  II  gives  the  similar  comparison  for  the  case  of  100%  overstrain. 

TABLE  I.  REDISTRIBUTION  OF  Oq  DUE  TO  TWO  ID  NOTCHES 
IN  A  CYLINDER  WITH  60%  OVERSTRAIN 


S  1  1  eg  in  Ksi  Corresponding  to  1 

1  1  1  (1  Ksi  -  6.9  MPa)  I 

T 

Node 

1  vu  r 

Thermal  Load 

mi 

Sr  +  ST| 

1 

(Fig.  4) 

1  1 

(1) 

[HI 

(4)  I 

T 

1  1 

1  r 

r 

1 

1 

1 

1  1.250  I 

-210.20 

1  -  38.91  I 

-171.28  | 

-210.191 

1 

2 

1  1.3067  | 

-146.75 

1  -  20.96  I 

-125.79  I 

-146.751 

1 

3 

1  1.3633  I 

-  95.87 

1  -  4.52  | 

-  91.41  | 

-  95.931 

1 

4 

1  1.420  I 

-  52.20 

1  10.63  | 

-  62.84  I 

-  52.201 

1 

5 

1  1.480  | 

-  14.37 

1  25.46  I 

-  39.80  I 

-  14.341 

1 

6 

1  1.540  I 

13.16 

1  39.20  | 

-  26.05  1 

13.151 

1 

7 

1  1.600  1 

34.10 

1  51.98  I 

-  17.84  1 

34.141 

1 

8 

1  1.7333  1 

45.85 

I  47.29  | 

-  1.46  1 

45.831 

1 

9 

1  1.8667  l 

60.78 

1  43.57  i 

17.23  1 

60.801 

1 

10 

1  2.0  | 

77.27 

1  40.57  | 

36.64  I 

77.211 

1 

88 

1  1.0  1 

-116.63 

1  -143.95  I 

26.92  1 

-117.031 

1 

89 

1  1.0833  1 

-  82.39 

I  -103.08  I 

20.98  1 

-  82.101 

1 

90 

1  1.1667  I 

-  53.89 

j  -  68.56  I 

15.03  1 

-  53.531 

1 

91 

1  1.250  I 

-  29.76 

1  -  38.91  i 

9.12  I 

-  29.791 

1 

92 

I  1.3067  | 

-  14.11 

1  -  20.96  I 

6.55  I 

-  14.411 

93 

1  1.3633  1 

-  0.80 

1  -  4.52  | 

3.78  1 

-  0.741 

1 

94 

1  1.42  1 

11.69 

i  10.63  j 

1.04  1 

11.671 

1 

95 

1  1.48  1 

24.24 

1  25.46  j 

-  1.45  1 

24.011 

1 

96 

1  1.54  1 

35.41 

1  39.20  | 

-  3.85  I 

35.351 

1 

97 

1  1.60  I 

46.01 

1  51.98  i 

-  5.97  I 

46.011 

1 

98 

1  1.7333  I 

36.95 

!  47.24  | 

-  10.48  1 

36.811 

1 

99 

j  1.8667  I 

29.32 

1  43.57  | 

-  14.32  1 

29.251 

100 

1  2.0  I 

23.21 

1  40.57  | 

-  17.46  I 

23.111 
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TABLE  It.  REDISTRIBUTION  OF  00  DUE  TO  TWO  ID  NOTCHES 
IN  A  CYLINDER  WITH  IOOX  OVERSTRAIN 


T 

1 - r 

oe  in  Ksl  Corresponding  to 

- r 

1 

1  1 

(1  Ksi  -  6.9  MPa) 

1 

T 

Node 

1  r/a  1 

Thermal  Load 

1  Sr  T 

ST 

Sr  +  STI 

1 (Fig.  4) 

1  1 

(1) 

1  (2)  I 

(3)  I 

(4)  1 

T 

1 

1 

1  1 

1  1.250  1 

-270.79 

1  r 

1  -  57.43  1 

1 

-213.39  I 

1 

-270.821 

1 

2 

1  1.3067  | 

-163.97 

1  -  30.70  | 

-133.26  1 

-163.961 

1 

3 

1  1.3633  I 

-  85.48 

1  -  7.14  | 

-  78.41  I 

-  85.551 

1 

4 

1  1.420  I 

-  31.11 

1  13.84  | 

-  44.89  1 

-  31.051 

1 

5 

1  1.480  1 

9.74 

!  32.72  I 

-  22.74  I 

9.78| 

1 

6 

I  1.540  1 

41.27 

1  49.85  I 

-  8.62  1 

41.231 

1 

7 

1  1.600  1 

65.80 

1  65.49  I 

0.31  I 

65.801 

1 

8 

I  1.7333  I 

92.89 

1  79.89  | 

13.02  1 

92.911 

1 

9 

1  1.8667  I 

122.06 

1  93.21  | 

28.85  f 

122.051 

1 

10 

1  2.0  | 

151.66 

1  105.59  | 

46.06  1 

151.651 

1 

88 

1  1.0  1 

-133.33 

1  -166.54  | 

32.68  1 

-133.861 

1 

89 

!  1.0833  1 

-  98.89 

1  -123.99  | 

25.42  1 

-  98.571 

1 

90 

1  1.1667  I 

-  70.45 

1  -  88.15  | 

18.16  1 

-  69.991 

1 

91 

1  1.250  1 

-  46.36 

1  -  57.43  | 

11.09  I 

-  46.341 

1 

92 

1  1.3067  | 

-  24.06 

1  -  30.70  | 

6.39  1 

-  24.311 

1 

93 

1  1.3633  1 

-  5.69 

1  -  7.14  | 

1.58  I 

-  5.561 

1 

94 

I  1.420  I 

11.03 

1  13.84  | 

-  2.86  I 

10.981 

1 

95 

I  1.480  I 

28.38 

1  32.72  I 

-  6.60  I 

26.121 

1 

96 

1  1.540  I 

39.69 

1  49.85  | 

-  10.14  I 

39.711 

1 

97 

1  1.600  1 

52.41 

1  65.49  I 

-  13.22  I 

52.271 

1 

98 

1  1.7333  1 

64.00 

1  79.89  | 

-  16.21  1 

63.681 

1 

99 

1  1.8667  | 

74.47 

1  93.21  | 

-  18.87  | 

74.341 

1 

100 

1  2.0  | 

84.30 

1  105.59  | 

-  21.25  1 

84.341 

CRACKS  INITIATED  FROM  NOTCHES 

In  the  previous  example  of  two  notches  In  an  autofrettaged  tube,  two 
radial  cracks  are  Initiated  from  the  notches,  Figure  5.  The  stress  intensity 
factor  at  the  crack  tip  can  be  computed  by  the  method  of  thermal  simulation 
without  additional  complications.  The  same  finite  element  idealization  shown 
In  Figure  4(b)  Is  used  here.  However,  elements  (1)  and  (2)  are  enriched  crack 
tip  elements.  The  stress  Intensity  factor  for  this  problem  Is  «  -100.87 
Ksi/Tn  for  a  thermal  load  corresponding  to  60%  overstrain. 
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To  use  the  method  of  superposlt ion,  it  Is  logical  to  obtain  the  stress 
distribution  along  the  crack  line  (nodes  1  to  4  In  Figure  4)  in  the  notched 
autofrettaged  cylinder  before  the  crack  is  actually  introduced.  This  computa¬ 
tion  has  been  accomplished  In  the  previous  section.  The  results  for  a  60% 
overstrained  cylinder  are  shown  In  Table  l.  Along  the  edge  of  element  1, 
where  the  radial  crack  will  appear,  the  tangential  stress  components  oq  are 
-210.19  Ksi,  -146.75  Ksi,  -95.93  Ksi  and  -52.20  Ksi  at  nodes  1,  2,  3,  and  4, 
respectively.  Applying  this  distributed  crack  face  load  as  the  only  loading 
condition  to  the  finite  element  idealization  shown  in  Figure  4(b)  again,  with 
elements  (1)  and  (2)  being  enriched  crack  tip  elements,  the  APES  program  gives 
the  stress  Intensity  factor  Kj  -  -103.41  Ksi/Tn.  The  percentage  of  error 
between  this  result  and  that  obtained  from  thermal  loading  is  about  three 
percent.  The  higher  percentage  of  error  is  due  to  the  accumulation  of  errors 
In  the  repeated  finite  element  computations  in  the  method  of  superposition. 

In  the  last  example,  the  cracks  appeared  after  the  presence  of  the 
notches.  The  order  of  presence  was  maintained  in  the  previous  computation. 
However,  the  order  is  indifferent  based  on  the  principle  of  superposition  in 
the  linear  theory  of  elasticity.  If  both  the  notches  and  the  cracks  are 
assumed  to  be  present  at  the  same  time,  it  provides  an  alternative  way  for  the 
computation  of  the  stress  intensity  factor.  Using  the  same  finite  element 
mesh  shown  in  Figure  4(b),  the  stress  intensity  factor  is  computed  for  distri¬ 
buted  loads  applied  on  an,  nq,  and  the  edge  from  node  1  to  node  4  of  element 
(1).  The  loads  equal  the  normal  stresses  on  mn,  nq  and  on  the  edge  of  element 
(1)  which  we  obtained  from  equations  (3)  and  (1)  for  the  autofrettaged  cylin¬ 
der  without  geometrical  discontinuities.  The  stress  intensity  factor  thus 
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obtained  from  APES  is  Kj  •  -100.81  Ksl/Tn.  It  Is  In  excellent  agreement  with 
Kj  «  -100.87  Ksi/Tn  for  the  simulated  thermal  load. 

With  the  method  of  superposition,  it  Is  necessary  to  apply  the  distrib¬ 
uted  load  on  the  crack  face  which  causes  singular  stress  at  and  near  the  crack 
tip  using  the  APES  program.  The  thermal  simulation  method  on  the  other  hand 
gives  quite  accurate  stress  near  the  crack  tip.  For  the  combination  of  cracks 
and  notches,  the  tangential  thermal  stress  at  node  5,  Figure  4(b),  is  oq  * 

-215  Ksl.  The  hoop  stress  estimated  from  the  stress  intensity  factor  based  on 
equation  (10)  at  the  same  node  is  oq  »  -232  Ksl. 

CONCLUSION 

The  thermal  stress  simulation  of  residual  stresses  can  be  used  for  the 
computation  of  stress  Intensity  factors  for  cracks  and  of  stress  redistribu¬ 
tion  due  to  the  presence  of  cracks  and  notches  In  a  residual  stress  field. 
Numerical  results  have  shown  the  excellent  agreement  between  thermal  stress 
and  residual  stress  and  between  the  redistributions  of  thermal  and  residual 
stresses  caused  by  notches  and  cracks.  The  method  Is  simple  and  accurate  and 
has  some  advantages  over  the  classical  method  of  superposition  which  results 
In  singular  stresses  near  the  crack  tip  in  the  finite  element  computation. 
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